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We propose two nonlinear random walk models which are suitable for the analysis of both chemo- 
taxis and anomalous transport. We derive the balance equations for the population density for the 
case when the transition rate for a random walk depends on residence time, chemotactic substance 
and population density. We introduce the anomalous chemotactic sensitivity and find anomalous 
aggregation phenomenon. So we suggest a new explanation of the well-known effect of chemotac- 
tic collapse. We develop a non-Markovian "velocity-jump" model and obtain the superdiffusive 
behavior of bacteria with power law "run" time. 

PACS numbers: 



Continuous time random walks (CTRW) have been 
widely used in many fields including physics, chemistry 
and life sciences (see, for example, excellent reviews 
[3, 0])- Many biological and physical transport processes 
exhibit anomalous behavior for which walker mean- 
squared displacement increases as a fractional power p 
of time: < x 2 (t) >^ (subdiffusion: \i < 1; superdif- 
fusion: fj, > 1). The chemotaxis is a directed migration 
of cells population toward a more favorable environment. 
The microscopic theory of the movement of cells or or- 
ganisms is also based on the random walk theory (see, 
for example, [3 6]). Although chemotaxis has a long 
history and has been studied by researchers for many 
decades, there is a lack of literature on connection be- 
tween the anomalous random walk and chemotaxis the- 
ory. We should mention recent exception |8| where biased 
CTRW has been analyzed. One of the reasons for this 
gap is that the chemotaxis is essentially nonhomogeneous 
in space and time random process, while the standard 
anomalous CTRW model involves the spatial and tempo- 
ral invariance [l], Q • The main purpose of this paper is to 
set up the random walk models for both chemotaxis and 
anomalous transport. These models can be also used in 
other biological applications of CTRW and chemotaxis as 
anomalous search strategy [7[ , cancer cell dichotomy |9j , 
subdiffusion in spiny dendrites (lpj . embryogenesis and 
wound healing. Much of the recent literature on chemo- 
taxis has been concerned with movement of bacteria E. 
coli involving the runs and tumbles. The standard as- 
sumption in modeling is that "run" and "tumble" time 
intervals are exponentially distributed However it 
has been found experimentally [ill ] that the distribution 
of "run" time intervals might deviate significantly from 
exponential approximation. It has a power law and leads 
to superdiffusive behavior of bacteria. One of the pur- 
poses of this paper is to examine the effect of long-time 
tails on bacteria movement in terms of "velocity-jump" 
model. 

We start with "space-jump" random walk model in one 
space dimension. The cell performs a random walk as fol- 
lows: it waits for a random time at each point in space 
before making a jump to another point. The most im- 
portant characteristic of this walk is the transition rate 7 
for jumps at point x. The standard assumption in CTRW 



theory is that 7 depends on the residence time (age) r. 
This is a time interval between two successive jumps of 
the cell. The corresponding waiting time density (f>(t) is 

related to 7(7-) as 4>(t) = j(t) exp (^f* j(T)drj [l2j]- In 

chemotaxis theory the jump of cells occurs in response 
to a chemical signal [4]. Therefore the transition rate 7 
should depend on chemotactic substance (external sig- 
nal) S(x, t) and its spatial and temporal gradient. It also 
depends on macroscopic population density p{x,t). This 
dependence describes the coupling of the cells density 
and chemotactic substance and crowding effects. Thus 



l(r\x,t) = j(T\S(x,t), S(x,t), p(x,t),t). 



(1) 



We introduce the cell density £,(x,t,T) at position x at 
time t with the residence time t. The main reason for 
introduction of the structured density £ is to make a ran- 
dom walk Markovian. This idea has been used in [H, Il2r - 
[l5j . The density £ obeys the balance equation 



§ + |^= -(r|.V(,-./). <V(.r./). ,,(.,-./). (2) 



We use the initial condition £(x, 0,r) = p (x)S(t) for 
which the residence time of all cells at t = equals to 
0; po(x) is the initial density of cells. It is clear that the 
residence time t varies from to t. The condition at 
t = can be written as 



£(x,t,0) = / / "f{r\x, t)^(x — z, t, t)w(z\x — z, t)dTdz. 
Jr Jo 

(3) 

Here w(z\x,t) is the dispersal kernel for jumps z which 
also depends on chemotactic substance and its gradient, 
density p(x, t) and t 



w{z\x, t) = w{z\S{x,t), S{x,t), p(x, t),t). 



(4) 



It is assumed that w is independent from t. On the left 
hand side of ((3J we have a density of cells just arriving 
at point x at time t (zero residence time). On the right 
hand side of ([3]) we have an integration of the rate at 
which the cells with different age r arriving at position 
x at time t from the different points x — z. Our purpose 
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now is to derive the Master equation for the cell density 



p(x,t) = / £{x,t,r)dT (5) 
Jo 

Using the method of characteristics, we find from ([2]) that 

£(x, t, r) = £(x, t - r, 0)e- A-r 7(«-(t-r)k,«)d. i (6) 

Let us denote the density of cells just arriving at point x 
at time t by j(x, t) = £(x, t, 0).We substitute (JB]) into §S§ 
and take into account the initial condition for £. We get 

j(x,t)= / i (x — z, t) w (z\x — z, t) dz, (7) 
jr 

where the i(x,t) is the density of cells leaving the point 
x exactly at time t : 



i(x, t) = / j(x, u)(f>(x, t, u)du + po(x)<fi(x, t, 0), (8) 
Jo 



t u ) = = l(t _ u | t y- H j{ s -u\x, s )d s 

at 

(9) 

and ^(x, t,u) is the probability that a cell is trapped at 
point x from time u to t without executing a jump 



^(a;, t, u) = e 



— p~ fu 7( s — u\x,s)ds 



(10) 



This is an extension of standard survival function for 
a nonlinear and nonhomogeneous case when \I> depends 
on chemotactic substance S(x,t) and population density 
p(x, t). The balance equation for p (x, t) can be found by 
substitution of ((6|) into ((5]) 



p(x,t) 



j (x.u) $>(x,t, u)du + p Q (x) %(x,t,Q). (11) 



The system of balance equations 0, © and pT|) is a 
nonlinear generalization of classical CTRW renewal equa- 
tions [l|, and CTRW models for inhomogeneous and 
nonlinear media [TH, [13] ■ These equations can serve 
as a starting point for the analysis of both chemotaxis 
and anomalous transport for "space-jump" random walk 
model. If we differentiate p (x, t) in (fTT|) with respect to 
time, we obtain the nonlinear Master equation 



dp 
dt 



i (x — z,t)w (z\x — z,t) dz — i (x. t) . (12) 



Now we are in a position to analyze the chemotaxis and 
anomalous effects in more detail. First we consider the 
case when a cell performs a random walk in a stationary 
environment with the distribution of chemotactic sub- 
stance S(x). In this case j(t\x, t) = 71 (r| £(#)). The sur- 
vival probability ^ in (|10[) must be a function of r = t — u 
and can be written as 



(13) 



Using the Laplace transform in (J7J, ([8]) and (jlll) . we ob- 
tain 

i(x,t) = f K x (t-T)p( X ,T)dT, (14) 
Jo 

where K x (t) is the memory kernel defined by its Laplace 
transform 



K x {s) 



Hs\S(x)) 

i>( s \s(x)y 



(15) 



where s is the Laplace variable. Substitution of (|14p into 
(IT2l gives a generalized Master equation dp/dt = L 1 p 
with the operator L l : 



L xP = f f 

Jo Jr 



Kx-z (t — t) p (x — z, r) to (z|x — z, i) dzdr 



K x (t — r)p (x, t) dr. 



(16) 



The case when the dispersal kernel w(z\x,t) depends on 
chemotactic substance S has been considered by Lang- 
lands and Henry [8j. It has been pointed out by Er- 
ban and Othmer that movement of bacteria in favorable 
environment is determined by chemokinesis rather than 
chemotaxis. In most cases the bacteria or cell "does not 
feel" a macroscopic gradient of S 6]. That is why it is 
more important to study the dependence of transition 
probability 7 on chemotactic substance S . To illustrate 
the general theory we use only a symmetrical dispersal 
kernel w(z) as a function of z. 

In a Markovian case, when 7 does not depend on the 
residence time variable r, we have ^(x,t) — e -71 ^™* 
and K x (s) — ji(S(x)). Under the diffusion approxima- 
tion, the Master equation (|16p takes the form 

ft =T& Msix))p(x > t)) > (17) 

where a 1 = J R z 2 w(z)dz. It is well known [3] that this 
equation can be rewritten as dp/dt + dJ/dx = with the 
flux of cells 



J 



dS _ o- a 7i(g(a)) dp 
2 dx' 



ox 



x(S(x)) 



a 2 dji 

Tits' 

where x is the chemotactic sensitivity. When the deriva- 
tive dji/dS is negative, the advection (taxis) is in the 
direction of increase in chemotactic substance. In gen- 
eral it follows from (TT6l) that cells flux is not local in 
time 



J 2 dx 



1 8K X (t - r) 



dS 



p (x, t) dr 



<j 2 r* 



dp (x,t) 
T) -^x~ dT - 



(19) 



Instead of \ we have a chemotaxis memory kernel 
dK x (t) I dS. Note that the memory kernel for the chemo- 
taxis flux is different form the memory kernel for diffusion 
term (compare to Q). 
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Let us consider the anomalous case when the waiting 
time PDF is heavy-tailed, such that the corresponding 
mean time is infinite. We assume that the longer cell 
survives at point x, the smaller the transition probabil- 
ity from x becomes. The rate j(t) is a monotonically 
decreasing function of residence time r. For example, if 
7 ( T ) = fi(S(x))/(P + t), it follows from JT3) that the 
survival function has a power-law dependence 



*(t\S(x)) = 



where /? is constant. Anomalous case corresponds 
to (jl(S(x)) < 1 | 0, when K x (s\S{x)) = 



'0 



tq is a parameter with units of time. 



The anomalous cell flux is 

a 2 dS d 1 



2 dx dS M s ^)) 



V 



i-MS(x)) 



p(x,t) 



2r, 



M (S(x)) 




V 



(S(x)) dp (x,t) 
dx 



(20) 



where the Riemann-Liouville fractional derivative 



V 



is defined p], HH as 



V 



l- M (S(z)) 



p(x,t) 



1 







p (x, u) du 



Y{n{S{x)))dt J (t - ' 

(21) 

It should be noted that the fractional time derivative of 
variable order p{x) has been considered in [16j |. When 
fi = const, we have a classical subdiffusion transport 
equation for which the mean squared displacement of cell 
increases with time as with fi < 1. 

Let us consider the aggregation phenomenon Q . In a 
Markovian case, in a finite domain with zero flux of cells 
on the boundary, there exists a stationary non-uniform 
solution of (fl~7)) [4| . The aggregation of cells is due to the 
fact that mean waiting time 7 ] ~ 1 (,5'(ir) is decreasing func- 
tion of the chemotactic substance S. In an anomalous 
case, the system is not ergodic and there is no stationary 
distribution. However, one can introduce the anomalous 
chemotactic sensitivity as a derivative of anomalous ex- 
ponent: = A 4 ' (S( x ))- When x M < 0, the cells will tend 
to aggregate where the exponent p is small. The anoma- 
lous flux (|2T)|) leads to p{x 7 t) — > S(x — x min ) as t — > oo. 
Here x m - m is the point where the anomalous exponent 
p (S(x)) has a minimum. It means that all cells aggregate 
into a tiny region of space forming high density system at 
the point x = x m in . This phenomenon can be referred to 
as anomalous aggregation. Similar results have been ob- 
tained in [18( for a simple two-state system. This effect 
is known in a literature as chemotactic collapse 0]. Here 
we suggest an explanation of this effect which is different 
from the classical one based on Keller-Segel equations. 
To prevent the occurrence of delta-distribution we need 
to take into account the crowding effect. In what fol- 
lows, we consider this effect by assuming that the tran- 
sition rate depends on both the residence time r and the 
population density p. 



If the transition rate 7 is independent of residence 
time r, then the system is Markovian. We assume 
that 7 depends on the time t and the density p (x, t) 
or non-stationary chemotactic substance S(x,t), that 
is j(r\x,t) — j2(p(x,t),t). Then we obtain i(x,t) — 
l2{p{x,t) 1 t)p{x,t). The nonlinear evolution equation for 
p is dp/ dt — L 2 p, where the operator L 2 is defined as 



L2P 



l2{p{x — z, t), t)p (x — z,t)w (z\x — z, t) dz 



-~/ 2 (p(x,t),t)p(x,t) . 



(22) 



Now let us consider the case when the transition proba- 
bility 7(r|x, f) depends both on the residence time r and 
the density p as follows 

1 (r\x,t)= 11 (r\S(x))+ l2 (p(x,t),t). (23) 

From (J7J, ([8} and (fTTj) . after lengthy calculations, we 
obtain 



i (x, t) 



K x (t - t) e~ £ i2{ P {x,s), s )ds p ^ T ) At 



+72 (p(x, t),t). 



(24) 



It turns out that the nonlocal term in (f2"4"]l involves the 
exponential factor with j2(p(x,t),t). Although 71 and 
72 are separable (see (|23l) ). the corresponding terms in 
(1241) are not separable. This is a non-Markovian memory 
effect. The generalized Master equation is dp/dt = hp, 
where 

Lp = J J K x - z {t-T)p{x-z,T) (25) 

xe -s:~<2( P (x~z,s),s)ds xw ( z \ x _ Zjt ) dzdT _ 

K x (t-r)p {x, t) e~ £ T*(/>(*,»).«)*'d T + L 2 p. 



It follows from here that Lp ^ L 1 p + L 2 p despite the 
fact that 7 = 7i+72- Similar phenomenon related to 
chemical reactions has been discussed in [3 Ell • The 
exponential factor with 72 in (|25p prevents an anomalous 
aggregation effect in a long-time limit. 

Let us consider now 1-D non-Markovian "velocity- 
jump" model for bacteria movement. The purpose is to 
get the superdiffusive behavior [ll|. The bacteria moves 
to the right with the velocity v + and reverses the direc- 
tion with the rate 7+. When the bacteria moves to the 
left with the velocity the turning rate is 7_. In gen- 
eral, the turning rate depends on run time t, on chemo- 
tactic substance S and macroscopic population density 
p(x,t) : j±(r\x,t) = j ± (T\S(x,t),S(x,t),p(x,t),t). Let 
£+(x, t, t) be the density of bacteria moving with veloc- 
ity v + with run time r. The corresponding density of 
organisms moving with the velocity V- is £—(x,T,t). In- 
tegration of £±(x, t, t) over the run time variable r gives 
the mean densities p±(x,t) — J Q £±(x,t,r)dT. The sys- 
tem of equations for £±(x, t, t) suggested by Alt [|| are 



dt dx 



dr 



-7±(rMK±. (26) 
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Initial conditions are £±(x, 0, t) = p±(x)8(t), where 
p±(x) are the initial densities. Boundary conditions at 
r = : 



£±(x,t,0) = / j T (T\S(x,t),p(x,t),t)£ T (x,t,T)dT. 
Jo 

(27) 

By using method of characteristics we solve (|2l))) and 
from (|27p after lengthy manipulations we find the non- 
linear system of equations for p± {x,t) and j± (x,t) = 
£±OM,0): 



P± (x,t) 



j± (x T v±(t — u), u) i &±(x, t, u)du 



+P± (xTv±t)^ ± {x,t,0), 



(28) 



j±(x,t) = - I j T (x ±v T (t - u),u)y T (x,t,u)du 
Jo 

-p Q T (x±v T t)y T (x,t,0). (29) 
Here we introduce the generalized survival function 

and its full derivative &± = ± u±^J-^ = 

— 7± (t\S(x, t),S(x, t), p(x,t),t\ If we differentiate 

p± {x, t) with respect to time, we obtain the system of 
nonlinear equations 

' !>± ± v± = j± (x, t) - j T (x, t) . (31) 



at 



dx 



If the switching rates 7+ and 7_ are in- 
dependent of run time r, then j± (x, i) = 
r y T (T\S(x,t), S(x,t), p(x,t),t)p T (x,t) . This hyper- 
bolic model has been studied by Hillen et al [5] . 

Let us illustrate the general theory by considering the 
case when v± = v and the run time PDF ip±( T ) — 



7±(t) exp (J Q r r )±{r)dr) behaves like 



+±{T) ~ (^) 



1+M 



H < 1, t -> 00. (32) 



The mean waiting time < r± >= TTp±(r)dT is in- 
finite. The experimental evidence of a power-law dis- 
tribution like (f3"2"]) has been reported in [ll|. We as- 
sume that v + = V- = v and all bacteria run in a 
positive direction initially. First we find the Laplace 

transform of (x(t)) : (x(s)) = —i ^ rfp ^' s - ) j ^ , where 

/o(fc, s) = p+(k, s) + p- (k, s) is the Fourier-Laplace trans- 
form of the total bacteria density of particles p = p + +p_. 
In the limit s 0, we find from (|28[) and (|29l) very un- 
usual result that (x(s)) ~ vTqS~ 2+ ^. It means that the 
average position of bacteria is not zero as it should be in 
Markovian case! It fact (x(t)) ~ vr^t 1 "^ as t —> 00. The 
spreading is slower than a ballistic motion (x(t) — vt) and 
faster than diffusion for < p < 0.5 (superdiffusion). 

In summary, we introduce two nonlinear CTRW mod- 
els which are suitable for the analysis of both chemotaxis 
and anomalous transport. We consider the case when 
the transition rate for a random walk depends not only 
on residence time, but also on chemotactic substance, 
its derivative and macroscopic population density. We 
manage to derive the balance equations for the popula- 
tion density and corresponding nonlinear Master equa- 
tions. We introduce the concept of anomalous chemo- 
tactic sensitivity as a derivative of anomalous exponent 
with respect to chemotaxis substance. We find the effect 
of anomalous aggregation when all bacteria tend to aggre- 
gate at the point where power-law exponent has a min- 
imum. So we suggest a new explanation of chemotactic 
collapse which is different from the classical one based on 
Keller-Segel equations. Motivated by experiment on run 
and tumble chemotaxis ill], we set up non-Markovian 
"velocity-jump" model and obtain the superdiffusive be- 
havior of bacteria with power law "run" time. 
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